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Effects of disorder on two-dimensional 7*2 topological insulator are studied numerically by the 
transfer matrix method. Based on the scaling analysis, the phase diagram is derived for a model 
of HgTe quantum well as a function of disorder strength and magnitude of the energy gap. In 
the presence of s z non-conserving spin-orbit coupling, a finite metallic region is found that 
partitions the two topologically distinct insulating phases. As disorder increases, a narrow- 
gap topologically trivial insulator undergoes a series of transitions; first to metal, second to 
topological insulator, third to metal, and finally back to trivial insulator. We show that this 
multiple transition is a consequence of two disorder effects; renormalization of the band gap, and 
Anderson localization. The metallic region found in the scaling analysis corresponds roughly 
to the region of finite density of states at the Fermi level evaluated in the self-consistent Born 
approximation. 
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Topological insulator is a distinct type of insulator, 
showing a metallic surface state. During the last few 
years an increasing attention has been on its Z2 ver- 
sion 1-3 ) which realizes without an external magnetic field 
and preserves time-reversal symmetry. This is contrast- 
ing to a more prototypical integer quantum Hall insu- 
lator. Another twist motivating the present study is a 
recent proposal on the concept of so-called "topological 
Anderson insulatorh (TAI). 4 ~ 6 ) which extends the limit 
of topological distinction beyond that of band insulator. 

In two spatial dimensions (2D) disorder plays a funda- 
mental role in electronic transport. This paper deals with 
a 2D version of the Z2 topological insulator (Z2TI), a 
variant of the model proposed by Bernevig, Hughes, and 
Zhang; hereafter referred to as BUZ. 7 ' The BHZ model 
has been proposed to describe a Z2TI implemented in 
HgTc/CdTc superstructure. 8,9 ^ An important element 
added in this study to the BHZ model is a Rashba-type 
spin-orbit coupling (SOC), which breaks spin-axial sym- 
metry along z-axis inherent to the original version of 
BHZ model. This s z non-conserving SOC changes the 
topological invariant characterizing insulator from Z to 
Z 2 . 10 ' n ) 

Effects of disorder on topological insulators have been 
extensively studied. 4 ~ 6 ' 12 ~ 18 ) Directly relevant to this 
Letter are refs., 12 ~ 14 ) which have treated the s z non- 
conserving cases of Z2TI, and refs. 4 ~ 6 ) dealing with TAI. 
The basic viewpoint on TAI on which we are based is due 
to ref., 6 ) i.e., spin-independent potential disorder brings 
about renormalization of the gap, both its magnitude and 
sign, allowing for conversion of an ordinary band insula- 
tor to a topologically non-trivial one. On the other hand, 
a metallic behavior induced by disorder is also expected 
in the s z non-conserving system we study. Since the ex- 
isting studies on TAI dealt with the s z conserving case, 



there was no discussion about how such metallic con- 
duction combines with the renormalization of band gap 
and localization effect to determine the phase diagram 
of disordered BHZ model with s z non-conserving SOC. 
This paper takes full account of such Rashba-type spin- 
flip terms incorporated in the BHZ model and deduces 
a phase diagram of the disordered and disorder-induced 
Z2TI. We demonstrate that disordered BHZ model shows 
as many as four transitions for a given value of trivial 
band gap as a function of disorder. As disorder increases, 
a narrow-gap topologically trivial insulator first turns to 
a metal, second to a topological insulator, third again to 
a metal, and finally back to the trivial insulator. Occur- 
rence of such "multiple transition" intervened by a finite 
metallic region is a distinct property, providing us also 
with a useful point of view on the understanding of TAI. 

We consider a tight-binding version of the BHZ model 
on 2D square lattice. We add Rashba-type SOC in or- 
der to generalize to include the electric field perpendic- 
ular to surface, which is induced by structural inversion 
asymmetry of the quantum well or gate electrode. With- 
out disorder, the Hamiltonian becomes block diagonal 
in momentum space and reduces to the following 4x4 
matrices, 7 ' 9,19 ' 20 -* 



H(k) 



h(k) r(fc) 
F(fc)t h*(-k) 



(1) 



The first (last) rows/columns correspond to real spin f 
(J.), i.e., s 2 = ±1. The 2x2 matrix F(fe) represents 
Rashba SOC, the explicit form of which will be deter- 
mined by symmetry arguments. The diagonal blocks, 
h(k) and its Kramers' partner h*(—k), are of 2D Dirac 
form: e.g., h(k) = d(k) ■ cr, where 

d(k) = (Asmk x , Asmk y , A — 25(2 — cos k x — cos k y )). 
cr is another set of Pauli matrices (than the real spin 
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s = {s x , s y , s z )), representing an orbital pseudo-spin. 
The parameter A, which is associated with the band 
gap, can be tuned by varying the thickness of HgTe 
layer. Then nontrivial phase with Z2 topological invari- 
ant v = 1 realizes when < A/B < 4 and 4 < A/B < 8, 
whereas trivial phase with v — corresponds to A/B < 
and A/B > 8. 

In the presence of a random potential, the total Hamil- 
tonian is defined only in real space, 



H 



~r+bty C, r 



(2) 



where r = (J, J) represents a site on 2D square lat- 
tice, and a = (1,0) and b = (0,1) are primitive lat- 
tice vectors with the lattice constant set to unity. Ex- 
plicit form of the hopping matrices t x and t y are de- 
termined as follows. We take the Rashba-type interac- 
tion, H R « (p x s) z = p x Sy — p y s Xl and adapt it to the 
BHZ model. The two components of the orbital pseudo- 
spin consist of Tq orbital with the total angular momen- 
tum j = 1/2, and one of Tg orbitals with j = 3/2 and 
j z = ±3/2. The form of H R is constrained due to various 
symmetry requirements. First, notice that hopping ma- 
trix elements transform under tt rotation with respect to 
the y-axis as, 



a/2=(a,j z \H R \0,-j z ) = 
= (0,j z \H R \ - a, -j z 



-a,-j z \(-H R )\0,j z ) 



(3) 



Here, particle-hole symmetry is assumed, implying that 
the matrix elements considered are independent of j z . 
Combined with translational invariance, eq. (3) guaran- 
tees that a is real. Second, invariance under 7r/2 rotation 
with respect to the z-axis requires, on the other hand, 



(a,j z \H R \0, -j z 



±i(b,j z \H R \0,- 
-(-a,j z \H R \0, 



Jz 



(4) 



where a positive (negative) sign corresponds to j z = 
±1/2 (±3/2). Note that ir/2 rotation gives rise to a fac- 
tor, exp(ij z 7r) = ±t. By combining eqs. (3) and (4), hop- 
ping matrices t x and t y are determined as, 



A 

= Ba z - 1— a a 



a 



t y = Ba z 



.A .a 

i — <x„ — i—a z 
2 y 2 



(5) 
(6) 



Thus off-diagonal coupling in cq. (1) has been determined 
as, 



r(fe) = ia diag [sin k x — i sin k 



y , sin k x + i sin k y 

]• (7) 
±3/2 states are 



Terms connecting j z = ±1/2 and j z 
neglected, for simplicity. 

For the impurity potential, we take a site-diagonal 
form. The local energy e r at site r is given by 



(A - 4B)a z (8)1+ diag 



(8) 



where I is 2 x 2 identity matrix in spin space, and strength 
of non-magnetic disorder is specified by Wr^ and Wr P ^ 
for Tq- and Ts-type orbitals, respectively. 7 ^ It is assumed 



that each impurity potential obeys uniform distribution 
in the period: \-W/2,W/2\. 

We have performed extensive numerical study on 
the localization length Al using the transfer matrix 
method. 21 ) For a given set of parameters (A,W), the 
localization length is estimated by varying the sys- 
tem size L. If Xl/L decreases (increases) with increasing 
L, the system is judged to be insulating (metallic). Here- 
after, we concentrate on the case of E = 0, where the 
Fermi level is located in the middle of band gap. The 
other parameters are fixed as A = B = 1. We iden- 
tify the localization characteristics from scaling behavior 
of the localization length A^, and determine the phase 
boundaries between metallic and insulating regions. 

Figure 1 shows representative data for different size 
of the system in the presence of s z non-conserving SOC 
with a = 0.5. In the upper panel, Xl/L is plotted as a 
function of W with A = — 1 fixed. At each value of W, Xl 
is evaluated for systems of different circumference L. It is 
clear that there are two values of W where Xl /L does not 
depend on L. This indicates a critical point that sepa- 
rates metallic and insulating regimes. In contrast with or- 
dinary phase transitions, the critical point is meaningful 
only at zero temperature, as an example of the quantum 
phase transition. As W increases from zero, the system 
becomes metallic at W = W\ ~ 7 up to W2 ~ 9.5. When 
W is further increased, Xl/L shows again a monotonic 
decrease as a function of L. Namely we observe a reen- 
trant behavior: localized — > metallic — > localized, with 
increasing W . The emergence of metallic state by the s z 
non-conserving SOC has already been studied by ref.12 
for the Kane-Mele model of graphene, which belongs to 
a topological insulator in the clean limit. In our case, 
however, the value A = — 1 means that the system is 
the trivial insulator in the clean limit. Namely, in con- 
tradiction to naive intuition, the disorder first converts 
insulator to metal, and then back to insulator. 

The lower panel of Fig. 1 shows Xl/L as a function 
of A with fixed W = 6. At each value of A, Xl/L is 
evaluated for systems of different circumference L. Over- 
all behavior of Xl/L looks similar to that of the upper 
panel. Namely, as A increases, the system evolves as: lo- 
calized — > metallic — > localized. Note that the first insu- 
lating phase: A<Ai^— 0.8 is of trivial nature (v = 0), 
whereas the second one with A > A2 ~ —0.5) is Z2 
nontrivial (v = 1). These different kinds are separated 
by a finite metallic region: Ai < A < A2. The second 
insulating phase with Z2-number v = 1 starts already 



at a negative value: A c 



-0.5. This result demon- 



strates that a Z2 topological insulator can be induced 
by introducing non-magnetic disorder to a clean trivial 
insulator. In recent literature, similar disorder-induced 
Z2-nontrivial phase has been discussed, 4 ~ 6 ' but without 
the intervening metallic region. 

By repeating such analysis for different values of A 
and W, and identifying critical points, we obtain the 
phase diagram of the system at the ground state. In the 
clean limit, the nontrivial phase appears in the region 
< A < 4 and 4 < A < 8. Only the region with A < 4 
is shown in Fig. 2 since the phase diagram is reflection 
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Fig. 1. (color online) Scaling behavior of as a function of W 
(upper panel) and A (lower panel). In the upper panel, the sys- 
tem with A = — 1 is topologically trivial in the clean limit. The 
lower panel shows the behavior of "Z2 topological Anderson in- 
sulator" , 4 ~ 6 ) since the region A > —0.5 is made topologically 
nontrivial by the effect of W. 
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Fig. 2. (color online) Phase diagram of disordered BHZ model in 
the presence (a = 0.5, blue triangles) and absence (a = 0, red 
circles) of s z non-conserving SOC. Lines connecting the data 
points are guide to eyes. 



symmetric about A = 4. The particle-hole symmetry is 
responsible for this property. For comparison, we also 
show the s z conserving case with a = 0. In this case the 
system is always insulating except along the transition 
line indicated by (red) circles. 

Let us focus on the s z non-conserving (a ^ 0) SOC 
effects. The triangle symbols (blue) show the transition 
points between metallic and insulating phases. It is clear 
that the metallic region emerges in the vicinity of the 
transition line at a = 0. Consequently, the two topologi- 
cally distinct insulating phases are always separated by a 
metallic region of finite width. Furthermore, v = 1 phase 
is extended toward A < by finite disorder. Unfortu- 
nately, in the weak-disordered region below W ~ 4 with 
a = 0.5, we have been unable to determine the transition 
point from the data up to 64 sites due to strong finite- 
size effect. The transition can occur, on the other hand, 
only at A = in the clean limit, and the transition line 
should be continuously connected. Therefore, with de- 
creasing W, it is reasonable to assume that the metallic 
corridor persists and converges to the point A = O.As 
a result, increasing the strength of disorder W, multiple 
transition occurs in such a way as: 

(i) Z2-nontrivial — > metallic — > ^-trivial for positive A; 

(ii) Z2-trivial — ► metallic — > Z2-nontrivial — > metallic — > 
Z 2 -trivial for -0.5 < A < 0. 

These multiple transition are to be contrasted with the 
case a = where intervening metallic phases are absent. 

It is natural to ask the origin of the reentrant behav- 
ior with negative A. In order to answer the question, 
we now turn to the density of states p(E) of the system. 
Let us first study how the region of nonzero p(0) is corre- 
lated to the region of metallic conductance. We follow the 
previous study 6 -* for the s z conserving case (a = 0), to 
employ the self-consistent Born approximation (SCBA). 
Note that interference of electronic wave functions, which 
is crucial for the Anderson localization, is beyond the 
scope of the SCBA. However, the SCBA does describe 
disorder-induced rcnormalization of A, whose sign dis- 
tinguishes whether the system is topologically trivial or 
not. 

Effects of disorder are taken into account in the SCBA 
as the self-energy E(-E) of the Green function G(k,E), 
which is the 4x4 matrix in our case. We decompose the 
self-energy matrix as E = So + E z ct z . Note that the self 
energy is a scalar in the (real) spin space due to time- 
reversal symmetry. The renormalization of parameters 
occurs as 

A -> A(E) = A + Re E Z (E), (9) 

E^E(E) =E-Rc £(,(£), (10) 

where | A| represents the renormalized energy gap, and E 
is the shift of the energy. The SCBA gives the following 
self-consistency equation: 

= — 1 7— I (G(fc,i?)), (ii) 



12 



where (G(k, E)) = ((S-^(fc)-E( J B))" 1 ) is the disorder- 
averaged Green function. Because of particle-hole sym- 
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Fig. 3. (color online) Density of states p(0) (upper part) and 
renormalized mass A (lower part) as a function of W with fixed 
A = —0.2 for different values of a. In the SCBA, A's for various 
a vanish commonly at W ~ 2.89, which corresponds roughly to 
the peak of p(0). 



metry, Re Eo(0) and I m E z (0) vanish. Hence the Fermi 
level at E = is not shifted. The density of states is 
given by 

1 24 
p(E) = Im tr G{E) = Im Y, (E). (12) 

We solve Eq. (11) and derive Im E o (0) and Re E z (0). 
Near the metal- insulator transition, Im E z (0) is small 
and converges only slowly in the iterative method. There- 
fore we apply a root finding procedure called Steffensen's 
method only for Im Eo(0) in order to accelerate the con- 
vergence. 

Figure 3 shows the density of states p(0) (upper part) 
together with the renormalized gap A(0) = A + E z (0) 
(lower part) at the particle-hole symmetric point E = 0. 
Results at fixed A = —0.2 but with different values of a is 
given for comparison. In the upper part p(0) is plotted as 
a function of W , which shows in all cases a finite range of 
finite p(0). According to the SCBA, the density of states 
p(0) vanishes outside this range. The change of p(0) as a 
function of W with decreasing a seems continuous down 
to the limit a = 0. Both the width and the magnitude of 
the region with finite p(0) decreases with decreasing a, 
but remains finite at a = 0. 

The lower part in Fig. 3 shows that A(0) changes sign 
around W ~ 2.89, which corresponds roughly to the peak 
of p(0). Note that the bare value in the present case is 
A = —0.2, which corresponds to the trivial insulator with 
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Fig. 4. (color line) Density of states in the strongly disordered 
region for the same parameters as in Fig. 3. 



v = 0. Thus we find that the disorder drives the sys- 
tem into topologically nontrivial regime by renormaliz- 
ing A(0) to a positive value. The change of A(0) depends 
only slightly on a as shown in Fig. 3. Correspondingly, 
p(0) becomes finite in a finite range of W on both sides 
of the band inversion point (A = 0). 

As W is further increased, p(0) vanishes and stays null 
until the second onset of the topological transition as 
shown in Fig. 4. Finite value for p(0) resumes slightly 
before the entrance to the metallic phase shown in Fig. 2. 
We have found that, after the second onset, p(0) never 
decreases but continues to grow with increasing W. In 
other words, the end of the metallic phase around W = 
10 shown in Fig. 2 cannot be reproduced by the SCBA. 
Although the density of states is finite for large W, the 
wave function is actually Anderson-localized. The SCBA 
cannot describe this localization property. Hence, it will 
safely be said that the reentrant behavior of the phase 
transition is partly understood by the renormalization of 
A(0) and the onset of p(0) within the SCBA, especially 
in the range of small W. 

In conclusion, we have derived the phase diagram 
of disordered BHZ model in the presence of s z non- 
conserving SOC. As a crucial effect of s z non-conserving 
SOC, we have demonstrated that metallic state intrudes 
each transition. While TAI is realized as a distinct state 
in (W 1 i?)-spacc, 4 ' it is connected continuously to clean 
topological insulator in (A, iy)-space. Furthermore, the 
origin of reentrant behavior is investigated in terms of 
the SCBA. As in the case of a = studied earlier, 6 ) 
we have found for a / that the disorder drives origi- 
nally negative A to positive. This change of sign causes 
the topological transition induced by disorder. However, 
as an important difference from the previous study, we 
find that metallic state always accompanies the topolog- 
ical transitions in the phase diagram. Actual 2D Z2TI 
systems such as HgTe quantum wells have a sizable s z 
non-conserving SOC, which can be tuned continuously 
by gate voltage. We hope that our results encourage ex- 
perimental detection of (multiple) topological transitions 
with metallic intruders. 
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